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00 ■ Abstract 

o . 

' We study a formulation of euclidean general relativity in which the dynamical 

variables are given by a sequence of real numbers A^, representing the eigenvalues of 
the Dirac operator on the curved spacetime. These quantities are diffeomorphism- 
invariant functions of the metric and they form an infinite set of "physical observ- 
bX). ables" for general relativity. Recent work of Connes and Chamseddine suggests that 

they can be taken as natural variables for an invariant description of the dynam- 
ics of gravity. We compute the Poisson brackets of the XnS, and find that these 
can be expressed in terms of the propagator of the linearized Einstein equations 
5^ I and the energy- momentum of the eigenspinors. We show that the eigenspinors' 

energy-momentum is the Jacobian matrix of the change of coordinates from the 
metric to the A„'s. We study a variant of the Connes-Chamseddine spectral action 
which eliminates a disturbing large cosmological term. We analyze the correspond- 
ing equations of motion and find that these are solved if the energy momenta of 
the eigenspinors scale linearly with the mass. Surprisingly, this scaling law codes 
Einstein's equations. Finally we study the coupling to a physical fermion field. 
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1 Introduction 



In this paper, we present a novel formulation of euclidean general relativity. More pre- 
cisely, we study a theory that approximates general relativity at large scale while modi- 
fying it at short scale. The theory is characterized by the fact that it is not formulated 
strictly as a field theory. The dynamical variables are not fields: they are an (infinite) 
sequence of real numbers A„, n being an integer. The relation between the usual repre- 
sentation of the gravitational field by means of the metric tensor g^jj{x), and the repre- 
sentation in terms of the A„'s is given by the fact that the A„'s are the eigenvalues of the 
Dirac operator D defined on the spacetime geometry described by g^jj{x). 

Thus, the general idea is to describe spacetime geometry by giving the eigen-frequencies 
of the spinors that can live on that spacetime. This is in the spirit of the well known 
mathematical problem: "Can one hear the shape of a drum?" IQ. Namely the prob- 
lem of characterizing a two-dimensional shape by means of the spectrum of the laplacian 
defined on that shape. (Within some approximation, this spectrum gives the Fourier 
decomposition of the sound emitted by a drum with that shape.) 

The theory we study is implicitly contained in the recent work of Connes and Chamsed- 
dine . Alain Connes' exciting and ambitious attempt to unravel a microscopic noncom- 
mutative structure of spacetime 0, ^, ^, |[ has generated, among several others, also the 
idea that the Dirac operator D encodes the full information about the spacetime geometry 
in a way usable for describing gravitational dynamics. First of all the geometry can be 
reconstructed from D. More precisely it can be reconstructed from the (normed) algebra 
generated by D and by the smooth functions / on spacetime. If x and y are spacetime 
points, then their geodesic distance can be expressed in terms of as 

d{x,y) = supf{\f{x)-f{y)\ : ||[A/]||<1}. (1) 

The sup is over all functions / whose commutator with D has norm less than one (here / 
and D are viewed as operators acting on the Hilbert space of the spinors on the manifold, 
and the norm is the natural operator norm). Secondly, there is a natural way of giving 
the dynamics in terms of D: The Einstein-Hilbert action with a cosmological term is 
approximated by the trace of a very simple function oi D P], as we recall below. These 
results suggest that one can take the Dirac operator D as the object representing the 
dynamical field, and try to develop a dynamical theory for D. 

In fact, in p|, the powerful machinery of non- commutative geometry is used to ele- 
gantly encode the Yang-Mills structure of the standard model into a non- commutative 
component of the spacetime geometry. Accordingly, (a generalized operator) D codes the 
gravitational field and the Yang-Mills fields, plus, as a very remarkable bonus, the Higgs 
fields as well. However, even independently from these results on the standard model, 
we think that the idea of encoding the dynamical field into D is very interesting also for 
gravity alone. This is the idea we pursue here. We take the purely gravitational compo- 
nent of the Connes-Chamseddine theory only and we remain in the regime of conventional 
commutative geometry. 

More in detail, we consider the idea of encoding the gravitational field into D in its 
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simplest form: we take its eigenvalues A„ as the dynamical variables of general relativity. 
The main reason for which we think this is interesting is that these eigenvalues form an 
infinite family of diffeomorphism invariant observables. Objects of this kind have long 
been sought for describing the geometry. General relativity teaches us that fundamental 
physics is invariant under active diffeomorphisms. Physically, this means that there is 
no fixed nondynamical structure with respect to which location or motion could be de- 
fined ||. Starting with Peter Bergamann's pioneering work a fully diffeomorphism 
invariant description of the geometry, free from the gauge redundancy of the usual for- 
malism, has long been sought M, but without much success so far. Such a description 



would be precious, in particular, for quantum gravity |]I0|- Now, the A„'s are precisely dif- 
feomorphism invariant quantities. (In fact, the invariance is only under diffeomorphisms 
which preserve the spin structure; however, only large diffeomorphisms can change the 
spin structure.) 

As a first step towards using these ideas in classical and/or quantum theories, we de- 
rive an expression for the Poisson brackets of the Dirac eigenvalues. We obtain this result 
by using the covariant formulation of the phase space of general relativity described in 
[TT[| and by extending a technique developed in [1^, where the eigenvalues of the three- 
dimensional Weyl operator, invariant under spatial diffeomorphism only, where analyzed. 
Surprisingly, we find that the Poisson brackets of the eigenvalues can be expressed explic- 
itly in terms of the energy-momentum tensors of the corresponding Dirac eigenspinors. 
These tensors form the Jacobian matrix of the change of coordinates between metric and 
eigenvalues. The brackets are quadratic in these tensors, with a kernel given by the prop- 
agator of the linearized Einstein equations. The energy-momentum tensor of the Dirac 
eigenspinors provides therefore a key tool for analyzing the representation of spacetime 
geometry in terms of Dirac eigenvalues. 

We also study a variant of the Chamseddine-Connes spectral action 0. In its simplest 
version, this action is a bit unrealistic because of a huge cosmological term. This term 
is disturbing not only phenomenologically, but also because it implies that the small- 
curvature geometries for which the spectral action approximates the Einstein-Hilbert ac- 
tion are not solutions of the theory. We study a version of the spectral action which 
eliminates the cosmological term. We close by analyzing the equations of motion, derived 
from this action. These are solved if the energy momenta of the high mass eigenspinors 
scale linearly with the mass. This scaling requirement approximates the vacuum Ein- 
stein equations. Thus, we obtain a representation of Einstein equations as a scaling law. 
Finally, we briefly describe the coupling of a Dirac fermion. 

Our results suggest that the Chamseddine-Connes gravitational theory can be viewed 
as a manageable theory, possibly with powerful applications to classical and quantum 
gravity. The theory reproduces general relativity at low energies; it is formulated in 
terms of fully diffeomorphism invariant variables; and, of course, it prompts fascinating 
extensions of the very notion of geometry. 



A condensed version of some of the results presented here has appeared in |TB|. An 
extension of these ideas to supergravity has been considered in [jl4|. 
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2 General Relativity in terms of Eigenvalues 



Consider Euclidean general relativity on a compact 4d (spin-) manifold M without bound- 
ary. For definitness, let us assume that M has the topology M = S*^ x 5^ (the manifold 
that may be relevant in studying the finite temperature quantum properties of a compact 
Robertson- Walker universe). The metric field is g^v{x) = E^^{x)Eyi{x), with E^^{x) the 
tetrad fields. Indices /i = 1, . . . , 4 are curved while / = 1, . . . , 4 are internal euclidean, 
raised and lowered by the Kronecker metric The spin connection tu^j is defined by 
d[fj,El^ = Lo^^jE^^, where square brackets indicate anti-symmetrization. The dynamics is 
determined by the Einstein-Hilbert action 

^£;H[i?] = ^ Jd'x^R, (2) 

where g and R are the determinant and the Ricci scalar of the metric respectively, G is 
the Newton constant. We put the speed of light equal to one. The Planck constant does 
not enter our considerations, which are purely classical. 

In spite of a widespread contrary belief, the phase space is a spacetime covariant notion: 
it is the space of the solutions of the equations of motion modulo gauge transformations 
13]. Here, gauge transformations are 4d diffeomorphisms and local rotations of the tetrad 
field. Thus, the phase space T of general relativity is the space of the tetrad fields that 
solve Einstein equations, modulo internal rotations and diffeomorphisms. Equivalently, it 
is the space of the Ricci flat "4-geometries" . We denote the space of smooth tetrad fields 
as £, the space of the tetrad fields that solve the Einstein equations as S (Solutions) and 
the space of the orbits of the gauge group in £^ as ^ (4- Geometries). The phase space T 
is the space of the orbits of the diffeomorphism group in S, or, equivalently, the subspace 
of Q determined by the Einstein equations. By definition, observables are functions on F 

i- 

We now define an infinite family of such observables. Consider the Hilbert space Ti of 
spinor fields ip on M. The scalar product is 



(^^,0) = J d^x^ V^(x)0(x), (3) 

with bar indicating complex conjugation, and the scalar product in spinor space being the 
natural one in C^. With 7^ being the (Euclidean) hermitian Dirac matrices, the curved 
Dirac operator is 

D = tyE^{d, + iu,jK7'7'')- (4) 

The operator D is a self-adjoint operator on 7i admitting a complete set of real eigenvalues 
A„ and "eigenspinors" ipn- The manifold M being compact, the spectrum is discrete 

Dlpri = \n'4'n , (5) 

The eigenvalues are labeled so that A„ < A„+i, with repeated multiplicity. Here n is 
integer (positive and negative) and we choose Aq to be the positive eigenvalue closest 
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to zero. For simplicity we assume that there are no zero modes. The eigenvalues have 
dimension of an inverse length. 

Notice that the Dirac operator depends on the gravitational field E; so do, of course, 
its eigenvalues as well. We indicate explicitly this dependence by writing D[E] and A„[i?]. 
The latter defines a discrete family of real- valued functions on A„ : E i — > A„[i?]. 
Equivalently, we have a function A from S into the space of infinite sequences 

Ej^ix) ^ A„[E]. 

The image X{S) of S under this map is contained in the cone A„ < A„+i of R°°. The 
functions A„ are invariant under 4d diffeomorphisms and under internal rotations of the 
tetrads. Therefore they are gauge invariant and they are well defined functions on ^. In 
particular, they are well defined on the phase space of general relativity F: thus, they are 
observables of pure general relativity. 

Two metric fields with the same collection {A„} are called "isospectral" . Isometric E 
fields are isospectral, but the converse needs not be true [jl], 0. Therefore A might not be 



injective even if restricted to Q. The A„'s might fail to coordinatize Q: they might fail to 
coordinatize F as well, although whether or not this happens is not clear to us. However, 
they presumably "almost do it". Following Connes 0, it is tempting to consider the 
physical hypothesis that isospectral Ricci fiat 4-geometries, are physical indistinguishable 
("Spectral hypothesis"). 



3 The Poisson Brackets 

A simplectic structure on F can be constructed in covariant form ||TT|. A vector field X 
on S can be written as a differential operator 

.Y = /d^xA»|Bl ^ (7) 

where is any solution of the Einstein equations for the tetrad field, linearized 

over the background E. A vector field [X] on F is given by an equivalence class of such 
vector fields X, modulo linearized gauge transformations of X^(a;). A linearized gauge 
transformation is given in tetrad general relativity by 

E ^ E + 6ff^^E (8) 
5,,,^El{x) = C,E'^{x)+p\x)fJKE^{x), (9) 

where w is a vector field (generating an infinitesimal diffeomorphism) , p is an 5*0(4) Lie 
algebra valued scalar field (a = 1, ... ,6), generating infinitesimal 4d rotations, C is the 
Lie derivative, and faK the generators of the vector representation S'(9(4). Two linearized 
field X and Y (around E) are gauge equivalent if 

Y = X + 6^^^E, (10) 
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for some couple {v,p). 

The simplectic two-form Q of general relativity is given by 

m, y) = ^ I n, {Xl y/) e-u. e-P^'' (11) 

where (X/ F/) = (X^ F/ - X/). From now on we put 327rG = 1. Both 

sides of ([TTl) are functions of i?, namely scalar functions on S; this E is used to transform 
internal indices into spacetime indices. Here S : a i — > x{a) is chosen to be a (compact 
non-contractible) three-dimensional surface, such that, topologically, M = Ex S*^ (so that 
it gives the single non trivial 3-cycle of M), but otherwise arbitrary, and Up its normal 
one-form. 

fl is degenerate precisely in the gauge directions; 

n{X,Y)[E] = iff Y = 5^,p^E, (12) 

thus it defines a non-degenerate simplectic two form on the space of the orbits, namely 
on r. The coefficients of Q form can be written as 

nZ{x,y) = I d^a Up [5{x,x{a))^J{y,x{a))] e^u^e^"^-. (13) 



Because of the degeneracy, VL has no inverse on S. However, let us (arbitrarily) fix a gauge 
(choose a representative field E for any four geometry, and, consequently, choose a field X 
in any equivalence class [X]). On the space of the gauge fixed fields, VL is non degenerate 
and we can invert it. Let P^/(a;, |/) be the inverse of the simplectic form matrix on this 
subspace, namely 

/ A / d'z Plt{x, y) nr^y, z) F^{z) = J d'z 6{x, z) 5^ 5^ F^{z). (14) 

for all solutions of the linearized Einstein equations F, satisfying the gauge condition 
chosen. Integrating over the delta functions, and using (|13[), we have 

jj'a Up [Plt{x,x{a)PpF^{x{a))] 6^^„6— = ^/(x). (15) 

This equation, where F is any solution of the linearized equations, defines P, in the chosen 
gauge. Then, we can write the Poisson bracket between two functions f,g on S as 

(16) 

If the functions / and g are gauge invariant, namely well defined on F, then the r.h.s of 
(p!^) is independent of the gauge chosen. But equation (|1^) is precisely the definition of 
the propagator of the linearized Einstein equations over the background E, in the chosen 
gauge. For instance, let us choose the surface S as = and fix the gauge with 

Xt = l, X,^ = 0, Xi = l, X: = 0. (17) 
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where a = 1, 2, 3 and i = 1, 2, 3. Then equation (|T5]) becomes 



1 



;i8) 



1 ^2 g^j^j-^ j. _ propagator can be 



where we have used the notation x 
easily recognized. 

Next, we need the Jacobian of the transformation from metric to eigenvalues. The 
variation of A„ for a variation of E can be computed using standard time independent 
quantum mechanics perturbation theory. For a self-adjoint operator D{v) depending on 
a parameter v and whose eigenvalues A„(f ) are non-degenerate, we have 



d\n{v) 
dv 



(19) 



This equation is well known for its application in elementary quantum mechanics. It 
can be obtained by varying v in the eigenvalue equation for D{v), taking the scalar 
product with one of the eigenvectors, and noticing that the terms with the variation of 
the eigenvectors cancel. We now apply this equation to our situation, assuming generic 
metrics with non- degenerate eigenvalues (we refer to |jT^ for the general situation). We 
want to compute the variation of Xn[E] for a small variation of the tetrad field E. Let 
Ej^[x) be an arbitrarily chosen tetrad field and v a real parameter, and consider the one 
parameter family of tetrad fields E^ with components 

Ej^{x)^El{x)+vEl{x), (20) 

Under the standard definition of functional derivative, the variation 5\n[E\/ 5Ej^{x) of the 
eigenvalues under a variation of the tetrad, is the distribution defined by 



d'x E[{x) 



d\n[Ev\ 



6El{x) 



dv 



v=Q 



Using (0), we have 
Explicitely 



dv 



d\n{Ev] 

dv 



- = i^jjuiKW 1 \^n[Ev\). 

dv 

d%,/^]^4E,] ME,]. 



In w = we have 



dXn [E-v] 



dv 



v=0 



d^y Jg[E] i:4E] 



dD[E., 



dv 



(21) 

(22) 
(23) 

(24) 



v=0 



Assuming the tetrad fields are suitably well behaved, we can rewrite this equation as 

dXn [Ev] 
dv 



d 




v=0 ~ 


v=0 ^ 



d'y Jg[E,] ME] D[E,] ^E] 
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d_ 

dv 



v=0 



d 

dv 



v=0 



dv 



^„[^] D[E] i,^[E] 



v=0 



d'y Jg[E,] ^^[E] D[E.,] 



dJg[Ey] 



dv 



ME] X4E] 2Pn[E] 



v=0 



d yJg[Ey] (?/'„-D[-E„]^/'n - Ki^ni^n)- 



(25) 



The last formula gives the variation of the action of a spinor field with mass A„ under a 
variation of the metric, (computed for the n-th eigenspinor of D[E]). But the variation of 
the action under a variation of the tetrad is a well known quantity in general relativity: it 
provides the general definition of the energy momentum tensor (here in tetrad notation) 
Tj'lx). Indeed, the Dirac energy-momentum tensor (density) is defined in general by 



Dirac ; 



(26) 



where 5'Dirac = / ^/g {'ipDip — Xipip) is the Dirac action of a spinor with "mass" A. (Since 
we have not put the Planck constant in the Dirac action, A has dimensions of an inverse 
length, rather than a mass.) See for instance [l^, where the explicit form of this tensor 



is also given. By denoting the energy momentum tensor of the eigenspinor ipn as Tnj {x) , 
we have obtained, from 



D and dl) 

5El{x) 



T 



nl 



(27) 



This equation gives the variation of the eigenvalues A„ under a variation of the tetrad 
EI^{x), namely the Jacobian matrix of the map A. The matrix elements of this Jacobian 
are given by the energy momentum tensor of the Dirac eigenspinors. This fact suggests 
that we can study the map A locally in the space of the metrics, by studying the space 
of the eigenspinor's energy-momenta. As far we know, little is known on the topology of 
the space of solutions of Euclidean Einstein's equations on a compact manifold. A local 
analysis in F would of course miss information on disconnected components of F. 

In order to avoid a possible confusion, we remark that the quantities A„ are invariant 
under diffeomorphisms, not under arbitrary changes of the metric or the tetrad fiels: the 
left hand side of (|27|) does not vanish in general. Finally, notice that the above derivations 



would go through for several other operators, beside the Dirac operator. In |T9[ a formula 
similar to (^) has been derived for any second order elliptic selfadjoint operator. 

By combining (^,0) and (^) we obtain our main result (in this equation we restore 
physical units for completeness): 



{A„, A™} = 327rG Jd'xjd''y r[„i'(x) P^;!{x,y) T^fM 



(28) 



8 



which gives the Poisson bracket of two eigenvalues of the Dirac operator in terms of the 
energy-momentum tensor of the two corresponding eigenspinors and of the propagator 
of the hnearized Einstein equations. The right hand side does not depend on the gauge 
chosen for P. 

Finally, if the transformation A between the "coordinates" E^^[x) and the "coordinates" 
\n is locally invertible on the phase space F, we can write the simplectic form directly in 
terms of the A^'s as 

VL = ilmn dXn A c/A^, (29) 

where a sum over indices is understood, and where flmn is defined by 

^mn T/i{x) T„,j{y) = n'^j{x,y). (30) 

Indeed, let dEj^{x) be a (basis) one-form on F, namely the infinitesimal difference between 
two solutions of Einstein equations, namely a solution of the Einstein equations linearized 
over E. We have then 

n = jd'xj d'y nrAx, y) dElix) A dE;^iy) 

= Jd^xJ d'yQmn T„^(x) dE^ix) A T„:;(y) dE',{y) 

= flmn dXn A dX^- (31) 

An explicit evaluation of the matrix Qnm would be of great interest. 



4 Action and Field Equations 

We now turn to the gravitational spectral action 0. This action contains a cutoff pa- 
rameter Iq with units of a length, which determines the scale at which the gravitational 
theory defined departs from general relativity. We may assume that /q is the Planck length 
lo ~ lO^^^cm (although we make no reference to quantum phenomena in the present con- 
text). We use also mo = l/^o, which has the same dimension as D and A„. The action 
depends also on a dimensionless cutoff function %(«), which vanishes for large u. The 
spectral action is then defined as 

SG[D] = KTr[xillD')]. (32) 

K is a multiplicative constant to be chosen to recover the right dimensions of the action 
and the multiplicative overall factor in (^. 



The action (|3^ ) approximates the Einstein-Hilbert action with a large cosmological 
term for "slowly varying" metrics with small curvature (with respect to the scale /q). 
Indeed, the heat kernel expansion allows to write (see [|1^ for a different derivation). 



SG{D) = {lor^foK f ^dx + {lor'h^ I R^dx + ... . (33) 
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The functions /o and /2 are defined by /o = ^ x{u)udu and /a = ^ x{u)du, 
the integrals being of the order of unity for the choice of cutoff function made. The other 
terms in ( pSj) are of higher order in Iq. 

The expansion (^) shows that the action is dominated by the Einstein- Hilbert 
action with a Planck-scale cosmological term. The presence of this term is a problem for 
the physical interpretation of the theory because the solutions of the equations of motions 
have Planck-scale Ricci scalar, and therefore they are all out of the regime for which 
the approximation taken is valid! However, the cosmological term can be canceled by 
replacing the function x with x, 

X{u) = x{u) - e^xi^u) , (34) 
with e << 1. Indeed, one finds /o = , ^ = (1 — e)/2. The modified action becomes 

Sg{D) = ^ j n^dx + ... . (35) 



'0 



M 



We obtain the Einstein-Hilbert action (§) by fixing 



(36) 



If Iq is the Planck length y/hG, then k = where h is the Planck constant, up to 
terms of order e. Low curvature geometries, for which the expansion ( |33|) holds are now 
solutions of the theory. Thus we obtain a theory that genuinely approximates pure general 
relativity at scales large compared to Iq. 

Let us now consider the equations of motion derived from this action. Following our 
philosophy, we want to regard the A„'s as the gravitational variables. The action can 
easily be expressed in terms of these variables: 

Sg[X]=^Y: xil'X)- (37) 



However, we cannot obtain (approximate) Einstein equations by simply varying (37) with 
respect to the A„'s: we must minimize (^) on the surface \{£)-, not on the entire R°°. In 
other words, the A„'s are not independent variables: there are relations among them, and 
these relations among them code the complexity of general relativity. We shall comment 
on these relations at the end of the paper. We can still obtain the equations of motion 
by varying Sq with respect to the tetrad field: 

_ 5Sg _ ^ dSG_ SXn _ ^ dxill^l) rj. /./ X /ooN 
5El{x) ^ dK SEf^ix) ^ dK ^ ^ 

Defining f{u) =: -^^x^u), (|38D becomes 

E/(^X) A„T4(x) = 0. (39) 
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These are the Einstein equations in the Dirac eigenvalues formahsm. 

The simplest choice for the cutoff function x(^) is to take it to be smooth and mono- 
tonic on with 

^^^) = |o if. > 1 + 5 ^40) 

where S << 1. Namely xi"^) is the smoothed-out characteristic function of the [0, 1] 
interval. With this choice, the action ( p2D is essentially (namely up to corrections of order 
6) simply {k times) the number of eigenvalues A„ with absolute value smaller that mo! 
Then the function f{u) vanishes everywhere except on two narrow peaks. A negative one 
(width 26 and height 1/26) centered at one; and a positive one (width 26 /e and height 
e^/26) around the arbitrary large number 1/e =: s >> 1. The first of these peaks gets 
contributions from A„'s such that A„ ~ mo, namely from Planck scale eigenvalues. The 
second from ones such that A„ ~ smo. Equations ( P^ are solved if the contributions 
of the two peaks cancel. This happens if below the Planck scale the energy momentum 
tensor scales as 

An(mo)P(l) Tn{mo)^i{x) = Xn{smo)P{s) Tn{smofl{x), (41) 

Here, p(l) and p(s) are the densities of eigenvalues of IqD"^ at the two peaks and the index 
n{t) is defined by 

h^li^i) = t. (42) 

For large n the growth of the eigenvalues of the Dirac operator is given by the Weyl 
formula A„ ~ \/27ry~^/^n^/^, where V is the volume. Using this, one derives immediately 
the eigenvalue densities, and simple algebra yields 

Tntix) = A„ lo To^(x) . (43) 

for n >> n{mp), where Toj{x) = Tn(mo)i{x) is the energy momentum at the Planck 
scale. We have shown that the dynamical equations for the geometry are solved if below 
the Planck length the energy-momentum of the eigenspinors scales as the eigenspinor's 
mass. In other words, we have expressed the Einstein equations as a scaling requirement 
on the energy-momenta of the very-high-frequency Dirac eigenspinors. 

We add a few considerations that shed some light on this scaling requirement. Notice 
that is formed by a term linear in the derivatives of the spinor field and a term 
independent from these. The latter is a function of {ip, E, d^E), quadratic in ip. 

Tni = ^nT'^^M^n + ^nj E, dE]. (44) 

If we expand the last term around a point of the manifold with local coordinates x, 
covariance and dimensional analysis require that 

Sni = coXnEl + ciRl + 02 REl + o(^yy (45) 
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T ^ 



for some fixed expansfon coefficients cq, ci and C2. Here i?^ is the Ricci tensor. 
To be convinced that terms of this form do appear, consider the foUowing. 

Tr ^iYRi^lJlK + ... 
Ri + ... 



(46) 



For sufficiently high n, the eigenspinors are focally approximated by plane waves in local 
cartesian coordinates. For these functions, if we double the mass the frequency doubles as 
well: if Am = i A„, then d^ipm = t dfj^ipn- It follows that in general the energy momentum 
scales as 



T„i = t 



+ 



ci Ri + C2 RE]] +0 



Ar^ 



(47) 



For large A„ we can disregard the last term, and therefore (HSf) requires that the second 



square bracket vanishes. Taking the trace we have i? = 0, using which we conclude i?^ = 0, 
which are the vacuum Einstein equations. Thus, the equations of motion are solved if the 
scaling requirement on the high mass eigenspinors' energy momenta is satisfied, and this 
requirement, in turn, yields vacuum Einstein equations at low energy scale. 



5 Matter couplings 

The spinors ipn that appear in the previous sections do not represent physical fermions. 
They are mathematical quantities used to capture aspects of the pure gravitational field. 
In particular, there is no sense in which they act back on the geometry. In oder to 
describe the physical system formed by a (classical) fermion field, say with "mass" (inverse 
wavelenght) m, interacting with general relativity, namely an interacting Dirac-Einstein 
system, we have to introduce a (physical) spinor field il){x). The action that governs the 
dynamics of a fermion field and its interaction with the gravitational field is the Dirac 
action 

^Dirac[^, E]= ( iijDi; - m^^) d^x = (^[D - m]^) (48) 



Therefore the Dirac-Einstein system is governed by the total action 

S[D,^] = Sg[D] + Sj,i,,,[D,^] = nTr [x{llD^)\ + {i,[D - m]^). (49) 

The natural thing to do in the context of the present formalism is to expand ip in the 
basis formed by the ipn- Namely to write 
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and to describe gravity in terms of the A„'s and the fermion in terms of its components 
Qn- The action becomes 

S[Xn,an] =J2 [f^ xillK) + i^n - rn)\an\'^ . (51) 

n 

The equations of motion are 

^ [2kII fillXl) A„ + |a„|2] Tn^x) = 0, (52) 

n 

{Xn -m) an = 0. (53) 

Eq. ( p2D corresponds to the Einstein equations with a source and ( |53D is the Dirac 
equation on a curved spacetime. Notice that the latter is algebraic, and it can be solved 
immediately. In order for a solution to exist there should exist an h such that 

Xn = m. (54) 

The solution is 

Qn = for all n n, 

Qn = a an arbitrary constant. (55) 

This is not surprising: a fermion field of mass m on a geometry characterized by the 
Dirac eigenvalues {A„} is given precisely by the eigenspinor ipn with eigenvalue equal to 
m. Using the solution of the Dirac equation, (|52D becomes 



h^lt E /(^oA.) A„ T„^(a;) = SttG \a\^Tn'i{x) (56) 



where we have used the value of k (pG)). /^From the results of the previous section, we 
recognize the left hand side as the Einstein tensor; the right hand side is the energy 
momentum tensor of the fermion. 



In the presence of matter, the scaling law (^) is altered. Using again the Weyl formula, 
we obtain with simple algebra 

Tn^i{x) = Xn lo [ro^(x) +« |apTji^(x)], (57) 

where 

a = . (58) 

Equation ( p7D is the "scaling law" form of the Einstein equations, modified by the matter 
source term. 

The extension of the theory to other conventional matter couplings should not be 
difficult, but we do not pursue it here. 
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6 Summary and perspectives 



We have discussed the possibihty of describing gravity by means of the Dirac operator 
eigenvalues. This possibihty has been opened by the recent work of Connes and Chamsed- 
dine. We think that these new ideas might open a novel window over the physics of space- 
time and find applications in classical and quantum gravitation. The main obstacle for a 
full development of this approach is its natural euclidean character, due to the fact that 
on a non-compact lorentzian spacetime the Dirac operator will not have discrete spectrum 



(but see pO[ for 'lorentzian' attempts). However, the present formalism might still find a 
natural application in quantum or in thermal quantum physics. 

We have elucidated some aspects of the dynamical structure of the theory in the A„ 
variables by computing their Poisson algebra. This is given in equation (p8|). Perhaps 
a quantum theory could be constructed in a diff-invariant manner by studying represen- 



tations of this algebra, as suggested by Connes and Isham [^. At present, the Poisson 
algebra is not given in closed form, since the right hand side of equation (^) is not 
expressed in terms of the A„ themselves. This difficulty could be faced by expanding 
the energy momentum tensors in terms of the eigenspinors themselves, as suggested by 
Hawkins 



We have also studied the equations of motion of (a version) of the Chamseddine- 
Connes spectral action. This action defines a theory that approximates general relativity 
at large scale, where it could be used as a tool in classical gravity. It would be interesting 
to explore the modifications to general relativity that it yields at short scale. We have 
given the expression for a fermion coupling in this formalism. We have found a puzzling 
and intriguing way of expressing the Einstein equations as a scaling law for the energy 
momenta of the ultra-high-frequency eigenspinors. 

The striking feature of the formalism discussed here is that the theory is formulated 
in terms of diffeomorphism invariant quantities. The A„'s are a family of diffeomorphism 
invariant observables in euclidean general relativity, which is presumably complete or 
"almost complete" (it could fail to distinguish possible isospectral and not isometric ge- 
ometries). It should be possible, at least in principle, to represent "physical observations" 
in pure gravity as a function of the A„'s alone. Another remarkable aspect of the spectral 
action is that it introduces a physical cutoff and an elementary physical length with- 
out breaking diffeomorphism invariance. The spectral action cuts off all high frequency 
modes, but it does so in a diffeomorphic invariant manner without introducing back- 
ground structures. Since the number of the remaining modes is determined by the ratio 
of the spacetime volume to the Planck scale, one may expect that a theory of this sort 
could have infrared but not ultraviolet divergences in the quantum regime. The quantum 
theory based on the spectral action is therefore very much worth exploring, we think. 

The key open problem, in our view, is to better understand the map A given in (P) 
and its range; namely the constraints that a sequence of real numbers A„ must satisfy, 
if it represents the spectrum of the Dirac operator of some geometry. This problem can 
be addressed locally (in phase space) by studying the tangent map to A. We have show 
that this tangent map is given explicitly by the eigenspinor's energy-momenta. One could 
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begin to study A around simple geometries, such as a flat 4-torus. On a more general 
ground, the constraints on the A„'s are presumably the core of the formulation of the 
gravitational theory that we have begun to explore here. They should be contained in 
Connes' axioms for D in the axiomatic definition of a spectral triple 0. The equations in 
these axioms capture the notion of Riemannian manifold algebraically, and they should 
code the constraints satisfied by the A^. Finding the explicit connection between the 
formalism studied here and Connes axioms' equations would be of great interest. 



Acknowledgments 

We thank R DePietri, J Frohlich, D Kastler and especially A Connes, for suggestions 
and conversations. GL tanks all members of DAMTP and in particular Prof. M. Green 
for the kind hospitality in Cambridge. Work supported by the Italian MURST and by 
NSF grant PHY-95-15506. 



References 

[1 

[2 



[3] 
[4] 
[5 



M Kac Ann math Monthly 73 (1966) 1; C Gordon, DL Webb, S Wolpert, Bull Amer 
Math Soc 27 (1992) 134. 



A Connes, Comm Math Phys 182 (1996) 155-176; AH Chamseddine, A Connes, Phys 
Rev Lett 24 (1996) 4868; Comm Math Phys 186 (1997) 731. 

A Connes, Non commutative geometry^ Academic Press 1994; J Math Phys 36 (1995) 
6194; 

A Connes, J Lott, Nucl Phys Supp B18 (1990) 295; in Proceedings of the 1991 Cargese 
summer school, edited by J Frohlich et al. Plenum, New York 1992. 

D Kastler, Rev Math Phys 5 (1993) 477; Rev Math Phys 8 (1996) 103; D Kastler, T 
Schiicker, Rev Math Phys 8 (1996) 205; "The Standard Models a la Connes-Lott" , 
|hep-th/9412T85| ; B lochum, D Kastler, T Schiicker, J Math Phys 36 (1995) 6232; 



'Fuzzy mass relations in the standard model", |hep-th/9507150| ; "Riemannian and 



Non-commutative Geometry in Physics", |hep-th/9511011 . L Carminati, B lochum 



D Kastler, T Schiicker, "On Connes' new principle of general relativity. Can spinors 
hear the forces of spacetime?" , Piep-th/9612228| ; CP Martin, JM Gracia-Bondfa, JC 



Varilly, The Standard Model as a noncommutative geometry: the low mass regime, 
Piep-th/ 9605001] , Phys Rep, to appear. 

[6] G Landi, An Introduction to Noncommutative Spaces and their Geometries, Springer 
1997. 

[7] PC Bergmann, Rev Mod Phys 33 (1961) 510. 



15 



[8] CJ Isham, "Topological and Global Aspects of Quantum Theory", in Relativity, 
Groups and Topology: Les Houches, edited by BS DeWitt and R Stora (North- 
Holland, Amsterdam 1984). C Rovelli, Classical and Quantum Gravity 8 (1991) 297; 
8, 317 (1991); "Half way through the woods", in The Cosmos of Science, J Ear- 
man and JD Norton editors, (University of Pittsburgh Press and Universitats Verlag 
Konstanz), to appear; "Localization in Quantum Field Theory: Which aspects or 
Quantum Field Theory are compatible with what we know about spacetime?", in 
the proceedings of the 1996 Boston Colloquium for Philosophy of Science: Reflec- 
tions on the Foundations of Quantum Field theory, Boston Center for the Philosophy 
of Science, Boston, March 1996, to appear. 

[9] PAM Dirac, Proc Roy Soc A249 (1958) 333; Phys Rev Lett 2 (1959) 368. A Komar, 
Phys Rev 111 (1958) 1182; K Kuchar, J Math Phys 22 (1981) 2640; D27 (1983) 2277. 
CG Torre, Phys Rev D48 (1993) 2373; Phys Rev Lett 70 (1993) 3525. ET Newman, 
PC Bergmann, Rev Mod Phys 29 (1959) 443. ET Newman, C Rovelli, Phys Rev Lett 
69 (1992) 1300. V Husain Class and Quant Grav 10 (1993) L233. JN Goldberg, J 
Lewandowski and C Stornaiolo, Comm Math Phys 148 (1992) 377. 



[10] CJ Isham, "Structural Issues in Quantum Gravity", Imperial/TP/95-96/07, |gr^ 
I qc/9510063|; to appear in the proceedings of the GR14. 



[11] A Ashtekar, L Bombelli, O Reula, in Mechanics, Analysis and Geometry: 200 Years 
after Lagrange, edited by M Francaviglia, Elsvier 1991. 

[12] R De Pietri, C Rovelh, Class and Quantum Grav 12 (1995) 1279. 

[13] G Landi, C Rovelh, Phys Rev Lett 78 (1997) 3051. 

[14] IV Vancea, Phys Rev Lett 79 (1997) 3121; "Euclidean Supergravity in Terms of Dirac 
Eigenvalues", ^r-qc/9710132| . 

[15] JL Lagrange, Mem CI Sci Math Phys Inst France (1808) pg 1; JM Souriau, Structure 
des systemes dynamiques, Dunod Paris 1970, pg 132, 146; P Chernoff, J Mardsen, 
Infinite Dimensional Hamiltonian Systems, Springer Lecture Notes in Mathematics, 
Vol 452, Berlin 1974; BF Schutz, Geometrical Methods in Mathematical Physics, 
Cambridge University Press, Cambridge 1980; C Crnkovic, E Witten, in Newton's 
tercentenary volume, edited by SW Hawking and W Israel, Cambridge University 
Press 1987. 

[16] PB Gilkey, Invariance theory, the heat equation and the Atiyah-Singer index theorem, 
pg 330, CRC Press 1995. 

[17] JP Bourguignon, P Gauduchon, Comm Math Phys 144 (1992) 581. 

[18] S Deser, P vanNieuwenhuizen, Phys Rev DIO (1974) 411. 

[19] V Moretti, Phys Rev D56 (1997) 7797. 



16 



[20] M Kalau, J Geom Phys 18 (1996) 349; E Hawkins, Comm Math Phys 187 (1997) 
471. 

[21] C Connes and C Isham, private communication. 
[22] E Hawkins, private communication. 



17 



